function h = plot(surface,varargin)
if nargin == 1
    varargin = {'bo-','linewidth',1};
end
h = plot(surface.x,surface.h,varargin{:});
set(gca,'FontSize',18);
xlabel('x (site)');
ylabel('h (site)');